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ABSTRACT 

We report on optical and infrared observations of the anomalous X-ray pulsar (AXP) 
IE 1048.1—5937, made during its ongoing X-ray flare which started in 2007 March. We detected 
the source in the optical I and near-infrared K s bands in two ground-based observations and obtained 
deep flux upper limits from four observations, including one with the Spitzer Space Telescope at 4.5 
and 8.0 /im. The detections indicate that the source was approximately 1.3-1.6 magnitudes brighter 
than in 2003-2006, when it was at the tail of a previous similar X-ray flare. Similar related flux vari- 
ations have been seen in two other AXPs during their X-ray outbursts, suggesting common behavior 
for large X-ray flux variation events in AXPs. The Spitzer flux limits are sufficiently deep that we can 
exclude mid-infrared emission similar to that from the AXP 4U 0142+61, which has been interpreted 
as arising from a dust disk around the AXP. The optical/near-infrared emission from IE 1048.1—5937 
probably has a magnetospheric origin. The similarity in the flux spectra of IE 1048.1—5937 and 
4U 0142+61 challenges the dust disk model proposed for the latter. 

Subject headings: stars: neutron — pulsars: individual (IE 1048.1—5937) — infrared: stars — X-rays: 
stars 



1. INTRODUCTION 

Supported by extensive observational studies over the 
past few years, it is generally believed that anoma- 
lous X-ray pulsars (AXPs) are magnetars — young neu- 
tron stars possessing ultra-high ~10 14 G magnetic fields. 
While AXPs are primarily known as X-ray sources, 
exhibiting a variety o f variability behavior related to 
their magnet ar nature ([Woods fc Thomps on 2006; Kaspil 
2007), we now know that 5 of 10 identified AXPs — 
4U 0142+61, IE 1048.1-5937, 1RXS J170849-400910, 
XTE J1810-197, and IE 22 59+586— are bright at near- 
infrared (NIR) wavelengths (|Woods fc Thompson! 120061 : 
McGill AXP online catalog 4 ). Among them, the AXP 4U 
0142+61 is detected from optical to mid-infrared (MIR) 
wavelengths, indicating a spectral energy distribution 
(SED) that can be described by a two-component model: 
one a power-law spectrum over optical VRI and NIR 
J bands, presumably arising from the magnetosphere, 
and one thermal b lackbody-like over the 2.2-8 /im range 
(|Wang et alj |2006), arising from a debris disk. The dis- 
covery of such an optical/IR SED from the AXP was un- 
expected, an d the two-component model rem ains contro- 
versial (e.g.. iDurant fc van KerkwiikH2006bl ). Ideally, in 
order to understand the optical/IR emission mechanism 
for AXPs and in particular determine whether the IR 
emission from 4U 0142+61 is unusual or generic among 
AXPs, SEDs of other AXPs must be observed. However, 
most known AXPs are highly extincted in the optical 
range and often extremely faint in the IR. 

Among the known AXPs, IE 1048.1-5937 is pecu- 
liar in that it has exh i bited two long-t erm X-ray flares 
(jGavriil fc Kaspill2004t iTam et allfeOOSft which have not 
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been seen from other AXPs thus far. Monitored with 
the Rossi X-ray Timing Explorer (RXTE), the first flare 
was found to start in 2002 April and lasted approxi- 
mately two years. At the begin ning of the flare, the 
NIR counterpart was discovered (jWang fc Chakrabartyl 
2002), and in observations more than a year later, the 
counterpart was found to be approximately 2 m agnitudes 
dimmer in JK S (|Durant fc van Kerkw iik 2005), suggest- 
ing large flux variations associated with the flare. Such 
flaring events thus provide an opportunity to study op- 
tical/IR emission from an AXP other than 4U 0142+61. 
In addition, during their X-ray outbursts, the AXPs 
IE 2259+586 and XTE J1810-197 both exhibited re- 
lated NIR brightening. In the firs t source, its NIR and X- 
ray flux changes were correlated (jTam et al.ll2004ft while 
in the latter, its NIR flu x did not exactly follow the de- 
cline of the X-ray flux (|Rea et alj 12004 ICamilo et al.l 
12007ft . Optical/IR observations of IE 1048.1-5937 dur- 
ing another X-ray flare would allow us to verify if the 
source exhibits similar emission behavior. 

For thes e reasons, after a second flare started on 2007 
March 24 (|Dib et al.ll2007ft . we made multi- wavelength 
observations of IE 1048.1—5937 with large ground-based 
telescopes and the Spitzer Space Telescope. In this paper, 
we report on our optical/IR observations during this on- 
going X-ray flare. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Ground-Based Observations 

2.1.1. NIR Imaging 

On 2007 April 04, we observed the IE 1048.1-5937 
field in K 8 band using Persson's Auxiliary Nasmyth In- 
frared Camera (PANIC: iMartini et al.l2004ft on the 6.5-m 
Baade Magellan Telescope at Las Campanas Observatory 
in Chile. The detector is a Rockwell Hawaii 1024x1024 
HgCdTe array, having a field of view (FOV) of 2' x 2' and 
a pixel scale of 0.125"/pixel. The total on-source expo- 
sure time was 15 min. During the exposure, the telescope 
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was dithered in a 3x3 grid with offsets of 10" to obtain 
a measurement of the sky background. The observing 
conditions were good, with 0.5" seeing in K s . 

On 2007 May 08, we observed the target field in the 
same band with the same telescope and camera. The 
total on-source time was 18 min, and the same observing 
strategy was used. During the exposure, the conditions 
were not stable, with K s seeing rapidly changing between 
0.6" and 0.8". 

2.1.2. Optical Imaging 

We observed the target field in the optical bands us- 
ing the 8.2-m Very Large Telescope (VLT; ANTU) at 
the European Southern Observatory (ESO) and the 8- 
m Gemini South Telescope at the Gemini Observatory. 
Both telescopes are located in Chile. 

The VLT observation was made on 2007 May 07, the 
same night that the second Magellan NIR image was ob- 
tained. The instrument used was the FOcal Reducer and 
low d ispersion Spectrograph fFORS2; lAppenzeller et all 
119981 ). which consists of two 2kx4k MIT CCDs and has a 
FOV of 6.8' x 6.8'. The CCD detectors were 2x2 binned, 
having a pixel scale of 0.25"/pixel. We obtained 15 3-min 
/-band images of the field, resulting in a 45 min total ex- 
posure. The telescope was dithered in a 3x5 grid, with 
offsets of 5"x2.5". The conditions were excellent, with 
0.5" seeing. 

The Gemini imaging observations were made on 2007 
June 24 and July 15. The instr ument was Gemini Multi- 
Object Spectrograph (GMOS; EookHiD IQQI). The 
detector array of GMOS consists of three 2048x4608 
EEV CCDs. The pixel scale is 0.073"/pixel, while we 
used a detector binning of 2 pixels for the observations. 
In the first observation, we obtained 8 r' and 5 i' im- 
ages, with exposure times of 5 and 3 min, respectively. 
The telescope was dithered for the exposures, with off- 
sets of 10". The seeing was approximately 0.9" during 
the observations. In the latter observation, we obtained 
16 5-min r' and 4 3-min i' images, with the same observ- 
ing strategy. The seeing was approximately 0.8". 

2.1.3. Data Reduction 

We used the IRAF data analysis package for data 
reduction. The images were bias-subtracted and flat- 
fielded. In addition, because the Gemini GMOS detec- 
tors have significant fringing in i 1 band, a fringe frame 
provided by the Gemini Observatory was used for sub- 
traction of the fringes in our i images. From each 
set of dithered images in one observation, a sky image 
was made by filtering out stars. The sky image was 
subtracted from the set of images, and then the sky- 
subtracted images were shifted and combined into one 
final image of the target field. A summary of the images 
we obtained is given in Table [T] 

2.2. Spitzer 4-5/8.0 \xm Imaging 

We also observed IE 1048.1-5937 on 2007 August 
10 with the Spitzer Space Telescope. The imaging in- 
strument used w as the Infrared Array Camera (IRAC; 
iFazio et al.ll2004| ). It operates in four channels at 3.6, 4.5, 
5.8, and 8.0 /im, while two adjacent fields are simultane- 
ously imaged in pairs (3.6 and 5.8 /im; 4.5 and 8.0 /im). 
We observed our target in the 4.5 (bandwidth 1.0 /im) 



and 8.0 /im (bandwidth 2.9 /im) channels. The detectors 
at the short and long wavelengths are InSb and Si:As 
devices, respectively, with 256x256 pixels and a plate 
scale of l'.'2/pixel. The field of view (FOV) is 5f2x5.'2. 
The frame time was 100 s, with 96.8 and 93.6 s effective 
exposure time per frame for the 4.5 and 8.0 /im data, re- 
spectively. The total exposure times in each observation 
were 53.2 min at 4.5 /im and 51.5 min at 8.0 /im. 

The raw image data were processed through the IRAC 
data pipelines (version S16.1.0) at the Spitzer Science 
Center (SSC). In the Basic Calibrated Data (BCD) 
pipeline, standard imaging data reductions, such as re- 
moval of the electronic bias, dark sky subtraction, flat- 
fielding, and linearization, are performed and individual 
flux-calibrated BCD frames are produced. In the post- 
BCD (PBCD) pipeline, radiation hits in BCD images are 
detected and excluded, and BCD frames are then com- 
bined into final PBCD mosaics. The details of the data 
reduction in the pipelin es can be found in t he IRAC Data 
Handbook (version 3.0: lReach et al.ll2006h . 

3. RESULTS 

3.1. Ground-based Observations 

The counterpart to IE 1048.1-5937 was detected in 
both NIR observations made on 2007 April 04 and the op- 
tical observations made on 2007 May 07. We performed 
PSF-fitting photometry to measure the brightness of the 
source. The nearby field star X 5 was used for flux cali- 
brati on in the K s and / bands (jDurant fc van Kerkwijkl 
I2005l h The source's magnitudes were K s = 19.9 ± 0.1 
and I — 24.9 ± 0.2. We did not detect the source in the 
other observations. Using the X 5 star for flux calibra- 
tion (|Wang & Chakrabartvl 120021 ) . we derived 3cr limit- 
ing magnitudes for the Magellan K s and Gemini r' and 
i' images. The results are given in Table [TJ 

3.2. Spitzer /IRAC Observations 

We did not detect the source in the Spitzer IRAC 
4.5 and 8.0 /im images, and the derived 3cr flux up- 
per limits were 5.2 and 21.8 /iJy, respectively. The 
fluxes correspond to limiting magnitudes of 18.9 and 16.2 
mag, for the zero m agnitude fluxes of 179.7 and 64.1 Jy 
(jReach et al.ll2005ft at the Spitzer /IRAC 4.5 and 8.0 /im 
bands, respectively. These results are also given in Ta- 
ble [T] The source region has been observed previously 
with IRAC at the same wavel ength bands in 20 05, with 
no counterpart found either (|Wang et al.ll2007l ). Com- 
paring to the previous results, our 2007 upper limits are 
approximately 2 times deeper. 

4. DISCUSSION 

The optical/NIR counterpart of IE 1048.1-5937 was 
observed a few times during 2003-2006 after the dis- 
covery of the counterpart, and the so urce's brightness 
had been low (e.g., K s ~ 21.0-21.5; see lTam et aLll2008l 
for details). Compared to the flux measurements, par- 
ticularly those obtained in 2003 April and June when 
the counterpart was detected in both the optic al I and 
NIR JK S bands (jDurant &: van Kerkwijkl 120051 : also see 
Table [2]), our results clearly indicate an optical/NIR 
brightening during the 2007 X-ray flare. Indeed, as- 
suming Ay = 5.4 mag, the unabsorbed optical/NIR- 
to-X-ray flux ratios for the 2003 and our measurements 
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are nearly identical (see Figure [T]). In this comparison, 
the unabsorbed 2-10 keV X-ray fluxes at two epochs 
given in Table [5] are us ed, and Ay is e stimated from 
N H = 9.7 x 10 21 cm" 2 (|Tam et all [2001) by using the 



relati on A v = A H /1.79 x 10 21 cm" 2 (jPredehl k Schmittl 
1995). Also we note that even though our second K s 
measurement was an upper limit (K s > 20.1), it was 
only 0.1 mag smal ler than that obtai ned from a detection 
two days later by llsrael et al.l f2007) . Thus in these two 
sets of observations, the optical/NIR and X-ray fluxes 
appeared to be correlated. 
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Fig. 1. — Unabsorbed optical/IR flux and flux upper limits 
(diamonds) of IE 1048.1-5937 measured in 2007 April-August, 
normalized by the unabsorbed 2—10 keV X-ray flux obtained on 
2007 April 28 (see Table [2} . The dotted line indicates our two 
nearly simultaneous observations in the / and K s bands, and 
the dashed line indicates the opti cal and NIR measurements ob - 
tained in 2003 April and June (Durant & van Kcrkwiik 2005). 
For comparison, the unabsorbed optical /IR— to-X-ray flux ratios 
for 4U 0142+61 are also show n (squares; [Durant & van K erkwiikl 
2006b; Wang & Kaspi 2008), with the vertical and horizontal 
lengths of the bars indicating the flux variation range and ef- 
fective bandwidth at each b and, respectively. Ay = 3.5 
(Durant & van Kcrkwiik 2006a) is used for dereddening. The un- 
absorbed 2-10 keV X-ray flux used fo r 4U 01424-61 was 6.6X10- 11 
ergs s _1 cm - 2 (IGonzalez et al.|[2007l) . 



However, as shown bv lTam et al.l (|2008f ). the observed 
NIR brightening at the beginning of the 2002 X-ray flare 
preceded the pulsed X-ray flux peak by ~2 months. If 
we assume that the relation between the pulsed flux F p 
and total flux F u F p oc if- 54 ([Tarn et al.ll2008h . is gen- 
erally true for this AXP, the 2-10 keV unabsorbed to- 
tal flux derived at the NIR observation time would have 
been ~10 -11 ergs s _1 cm -2 , 3.6 times lower than the 
2007 peak flux. As a result, the 2002 NIR-to-X-ray flux 
ratios would be approximately 6 times higher than that 
shown in Figure[T](note K s — 19.4 in the 2002 detection). 
This would indicate that the NIR and X-ray fluxes were 
not correlated all the time. We note that since the F p 
and Ft relation is mainly based on the total flux mea- 
surements obtained during the 2007 flare and there were 
not such measurements for the 2002 flare, it is possi- 
ble that the total flux was much larger at the start of 
the 2002 flare, and not reflected in the RXTE pulsed 
flux measurements. In any case, we confirm that in this 
AXP, the large optical/NIR flux variations were related 
to the X-ray flares. Similar related NIR brightening has 



been seen in IE 2259+5 86 and XTE J1810-197 dur- 
ing th eir X-ray outbursts (|Tam et alj|2004t iCamilo et al.l 
suggesting common behavior for emission from 
AXPs. However for 4U 0142+61, while its X -ray flux has 
been relatively stable (|Gonzalez et al.l r2007f ) . its optical 
and NIR flux was found to hav e large, rapid variations 
(|Durant k van Kerk wiik 2006b). This limits the related 
NIR/X-ray behavior to only large X-ray flux variation 
events. 

It is instructive to compare IE 1048.1—5937 to 
4U 0142+61, since the latter has well measured fluxes 
from the optical to MIR. In Figure[TJ we compare the op- 
tical/IR fluxes and flux upper limits of IE 1048.1-5937 
to those of 4U 0142+61, with each set of data points 
normalized by the approximately contemporaneous un- 
absorbed 2-10 keV X-ray flux of each source. The optical 
VRI and NIR JHK S fluxes for 4U 0142+61 are from Du- 
rant & van Kerkwijk (2006b; because of relatively large 
flux variations from the source, we show the range of the 
flux found in each band in the figure) , the two MIR data 
points at 4.5 and 8.0 fim are the 2005 flux measurements 
in Wang & Kaspi (2008), and the unabsorbed 2-10 keV 
X-ray flux used is 6.6xl0~ 11 ergs s" 1 cm -2 , obtained 
on 2004 July 24 (jGonzalez et al.ll2007T) . As can be seen 
from the figure, while the flux ratios of IE 1048.1—5937 
are roughly two times larger than those of 4U 0142+61, 
the deep limit at 4.5 /jm from our Spitzer IRAC observa- 
tion is approximately 3 times lower than the detection of 
4U 0142+61. This strongly suggests that our target does 
not have similar MIR emission. This difference in the 
MIR is not very sensitive to uncertainties on the redden- 
ing to the two sources, since MIR emission is only weakly 
extincted by the interstellar med ium. The X-ray flu x 
from 4U 0142+61 has been stable (jGonzalez et al.ll2007h . 
A possible explanation could be the non-simultaneous 
X-ray flux of IE 1048.1—5937 we use for the compari- 
son. However, RXTE X-ray monitoring observations of 
the source have indicated tha t the pulsed flux ; which is 
found to trace the total flux (|Tam et al.ll2008D . has re- 
mained stable and high during our observations (Dib ct 
al. 2008, in preparation). In order to raise the 4.5 jum 
upper limit to the flux ratio of 4U 0142+61, the X-ray 
flux would have had to have been as low as llxl0~ 12 
ergs s _1 cm -2 , only 40% larger than the average quies- 
cent flux (7.7xl0~ 12 ergs s _1 cm~ 2 ) of the source and 
approximately on e third of those obtained in 2007 April 
(|Tam et al.ll2008l ). 

The lack of MIR emission similar to that of 
4U 0142+61 from IE 1048.1-5937 seems ro bust. In the 
dust disk model proposed for 4U 0142+61 (Wa ng et al.l 
2006), the MIR flux, arising from X-ray irradiation of 
the disk, is proportional to the X-ray flux of the pul- 
sar. Therefore, even though our deep 4.5 /xm limit was 
obtained during the X-ray flare, we can exclude the ex- 
istence of a similar disk around IE 1048.1—5937. How- 
ever, since a disk could be further away from the cen- 
tral source, deep observations at longer wavelengths are 
needed in order to exclude the existence of a disk more 
conclusively. The MIR non-detection may suggest that 
debris disks are not commonly found around AXPs and 
thus the putative disk in 4U 0142+61 is unique. Two 
other AXPs have also been observed by Spitzer with no 
detection; however, the derived upper limits are far above 
the MID-to-X-ray flux ratio of 4U 0142+61 (jWang et all 
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12001 . 

The optical and NIR emission from IE 1048.1-5937 
probably has a magnetosphcric origin, given that 
the flux spectrum is s imilar t o that of 4U 0142+61 
(I Wang fc Chakrabartvl 120021 : iDurant fe van Kerkwijk! 
I2005h . For the latter source, its optical emission is 
known to be pulsed at the spin period and has a 
pulsed fraction muc h higher than that in X-rays, ex- 
cludi ng a disk origin (|Kern fc Martiij l2002; Dhil lon et all 
I2005h . Details of how optical and NIR emission is 
produced in the magnetosphere of a magnetar are 
not known, though p ossible radiation mechanisms hav e 
been suggested (e.g., iBeloborodov fc Thompson! l2007f ). 
We note that in the disk model for 4U 0142+61 
(jWang et all [2006), the if -band flux primarily arises 
from the disk, not from the magnetosphere. There- 
fore, the simil arity in the if-band-to-X-ray flux ratios 
for the AXPs (|Durant fe van Kerkwi ik 2005), including 
IE 1048.1— 5937 as confirmed by our measurements, chal- 
lenges the disk model. The if -band emission thus should 
originate from the magnetosphere, which may be sup- 
ported by the fact that in the absence of X-ray vari- 
ability, strong if-band variability is seen in 4U 0142+61 



(jDurant fc van Kerkwiik1l2006bf) . 
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TABLE 1 

Optical/IR Observations of IE 1048.1-5937 



Observation Start 


MJD 


Telescope /Instrument 


Filter 


Exposure 


Seeing a 


Magnitude 1 " 


J/F„ c /10- 14 


Time (UTC) 








(min) 


(arcsec) 




(ergs s x cm ') 


2007-04-04 00:37 


54194.0 


Magellan/PANIC 


K s 


15 


0.54 


19.9±0.1 


1.8 


2007-05-07 23:24 


54228.0 


VLT/FORS2 


I 


15 


0.55 


24.9±0.2 


1.9 


2007-05-08 03:15 


54228.1 


Magellan/PANIC 


K s 


18 


0.66 


>20.1 


<1.5 


2007-06-24 01:18 


54275.1 


Gemini-S/GMOS 


r' 


40 


1.0 


>24.8 










%> 


15 


0.88 


>24.1 




2007-07-15 00:06 


54296.0 


Gemini-S/GMOS 


r' 


80 


0.92 


>25.6 


<6.6 








i' 


12 


0.74 


>24.5 


<5.4 


2007-08-10 10:04 


54322.4 


Spitzer/IRAC 


4.5 fim 


53.2 




> 18.9 


<0.42 








8.0 ^m 


51.5 




> 16.2 


<0.99 



a Full width at half maximum. The upper limits are 3er limiting magnitudes. Flux is dereddened with Av — 5.4, using the reddening laws for 
the optical and NIR bands given in Schlcgcl ct al. ( 1998) and for the Spitzer/IRAC data given in Indcbctouw ct al. (2003). 



TABLE 2 

Flux measurements of IE 1048.1—5937 used in the 

COMPARISON. 



Epoch 


Optical/NIR 


X-ray a Refs 




Magnitude 


vF v h 




2003 April 24 


J = 23.4 


6.6 


1 


2003 June 6 


/ = 26.2 


5.9 


1 


2003 June 7 


K s = 21.5 


1.9 


1 


2003 July 16 






11.5 2 


2007 April 28 






36.1 2 
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a X-ray flux is phase-averaged, unabsorbed, in the energy range of 
2-10 keV, and in units of 10 -12 ergs s -1 cm -2 .* 3 Dereddened flux 
(assuming Av — 5.4; in units of 10 — 15 ergs s" 1 cm -2 ) obtained 
at the tail of the first X-ray flare. 



